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ABSTRACT 

We calculate the colour averaged and colour singlet free energies of static quark anti- 
quark sources placed in a thermal gluonic heat bath. We discuss the renormalization 
of these free energies using the short distance properties of the zero temperature 
heavy quark potential. This leads to the definition of the renormalized Polyakov 
loop as an order parameter for the deconfmement phase transition of the SU(3) 
gauge theory which is well behaved in the continuum limit. 



1 Introduction 



Universal properties of the finite temperature phase transition in non-abelian SU(N) 
gauge theories generally are discussed in terms of the Polyakov loop, which is an or- 
der parameter for the confinement-deconfinement transition. In a somewhat loosely 
defined terminology the Polyakov loop, L, is said to be related to the free energy of 
a heavy, static quark put as a test charge into a thermal medium, L ~ exp(— F q /T). 
However, it is well known that the Polyakov loop calculated on the lattice is a priori 
not properly renormalized. Taking the continuum limit at fixed temperature will 
lead to a vanishing Polyakov loop expectation value even in the deconfined phase as 
there are divergent self-energy contributions to the free energy. A proper renormal- 
ization of the Polyakov loop thus is needed in order to relate it to the heavy quark 
free energy in the continuum limit which then can also be used to define a proper 
order parameter that survives the continuum limit and can, for instance, be used to 
construct effective actions for the confinement-deconfinement transition [[J . 

We suggest here that a proper renormalization of Polyakov loop can be ob- 
tained through the renormalization of the finite temperature free energy of a static 
quark anti-quark pair calculated at short distances. At short distances the quark 
anti-quark pair interacts through the exchange of gluons, which may be calculated 
perturbatively. More importantly we expect that this interaction is essentially tem- 
perature independent for separations r which are smaller than the average separation 
between partons in the thermal medium. At short distances the finite temperature 
heavy quark anti-quark free energy thus will be given by the zero temperature heavy 
quark potential. This allows to remove divergent self-energy contributions in the free 
energy, F qq (r,T), through a matching of its short distance behaviour to that of the 
heavy quark anti-quark potential, V qq (r), at zero temperature. After having done 
so also the large distance behaviour of F qq (r,T) is fixed. As this, in turn, is re- 
lated to the Polyakov loop expectation value, it allows to define a renormalized 
Polyakov loop, i.e. as usual no additional divergences will show up in the heavy 
quark free energy at finite temperature, once it has been renormalized properly at 
zero temperature. 

It should be obvious that the program outlined above equally well holds for the 
heavy quark free energy calculated in QCD, i.e. in the presence of light dynamical 
quarks. It will be of even greater importance in this case as F qq (r, T) will be finite 
for all temperatures and its large distance behaviour in the low temperature, chiral 
symmetry broken phase, will reflect the temperature dependence of the string break- 
ing energy. In fact, the conceptual approach we are going to present and analyze 
here has already been anticipated in the presentation of heavy quark anti-quark free 
energies given in Ref. f2|. 

This paper is organized as follows: In Section 2 we summarize basic definitions 
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for colour averaged and singlet free energies and outline our approach to renormalize 
the Polyakov loop expectation value. In Section 3 we present a detailed calculation 
of the heavy quark anti-quark free energy at short distances. These calculations 
are performed within the SU(3) gauge theory on lattices with temporal extent up 
to N T = 16 so that small distances in units of the temperature can be resolved, 
rT ~ 1/N T . We also point out the usefulness of the colour singlet free energy for 
the matching to the heavy quark potential at short distances and the determination 
of the asymptotic behaviour of free energies at large distances. In Section 4 we 
discuss the properties of the renormalized Polyakov loop. Section 5 contains our 
conclusions. 



2 Heavy quark anti-quark free energy 

The free energy of a pair of static quark anti-quark sources in a thermal medium 

Fqq{r T T) = -ln^fr^frLt)^) + c(T) , rT = \x - y\/N T , (1) 

is represented in terms of the Polyakov loop, which on an Euclidean lattice of size 
N% x N T is defined as a product of gauge field variables U X0) g, 



n 
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In Eq. |I] and in what follows we use the notation Tr = |Tr. Eq. [T] defines the 
colour averaged free energy up to an additive normalization constant, c(T), which 
is related to the self-energy of the quark and anti-quark sources. The Polyakov loop 
expectation value, (L) = (N~ 3 J2g TrL^), is closely related to the long (infinite) 
distance behaviour of the free energy, 

F oa (T)= lim F qq (r,T) = -T In | (L) \ 2 + c(T) T . (3) 

r— >oo 

It is obvious that we need to fix the normalization constant c(T) in order to give a 
physical meaning to this free energy. 

In the pure SU(3) gauge theory the expectation value of the Polyakov loop, 
provides an order parameter for the deconfinement phase transition; in the confined 
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phase exp(— Foo(T)/2T) vanishes, while it stays non-zero in the deconfined phase. 
We would like to interpret as the change in free energy due to the presence of two 
well separated quarks, each of which is screened by a cloud of gluons. In order to do 
so we have to clarify the normalization of F qq (r,T) which unambiguously fixes the 
T-dependent additive constant in Eq. [j]. Before doing so it, however, is worthwhile 
to discuss in a bit more detail the properties of F qq (r, T). 

In addition to the colour averaged potential defined through Eq. |I] we introduce 
the colour singlet potential ||, 

= - In ((Tr^4>) + c'(T) . (4) 

As it stands this correlation function is gauge dependent; for its evaluation we 
have to fix a gauge. However, it recently has been shown that a gauge independent 
definition of the singlet free energy can be given || and that this definition coincides 
with the definition given by Eq. |] in gauges which have a positive transfer matrix 
representation on the lattice. Calculating F\ from the Polyakov loop correlation 
function given in Eq. ^ in, for instance, the Coulomb gauge thus yields the singlet 
potential. 

The colour averaged as well as the singlet free energy will approach finite, non- 
zero values for all temperatures T > T c (and diverge below T c ). In fact, if the 
separation between the quark anti-quark sources gets large the relative orientation 
of the charges in colour space will not influence the screening of the individual 
charges. We thus choose the relative normalization of F± and F qq such that they are 
identical in the limit of large spatial separations. The colour averaged free energy 
defined in Eq. [TJ may then be represented by a thermal average over free energies of 
a qq-pair in singlet (Fi) and octet (F 8 ) representations, respectively |3|, ||, 

exp(-F $? (r,T)/T) = ~ex P (-F 1 (r,T)/T) + ~ exp(-F 8 (r, T)/T) . (5) 

At distances much shorter than the inverse temperature (rT << 1) the dominant 
scale is set by r, the running coupling will be controlled by this scale, g(r,T) ~ 
g(r, T — 0) = g(r), and will become small for r « 1/ 'Aqcd- In this limit the singlet 
and octet free energies are dominated by one-gluon exchange and are calculable 
within ordinary zero temperature perturbation theory, 

F 1 (r,T = 0) = -8F 8 (r,T = 0) + O(g i ) = -^(l + O(g 2 )) , (6) 

07IT \ J 

where F\{r,T = 0) = Vq q (r) is the static, zero temperature heavy quark potential. 
As the singlet potential is attractive and the octet potential is repulsive at short 
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distances it follows from Eq. [5] that in this limit the colour averaged free energy will 
be dominated by the singlet contribution. We then may deduce from Eq. |5| also the 
asymptotic short distance behaviour of F qq and Ti, 

]im(F qq {r,T) - Fi(r,T)) = Tin 9 for all T. (7) 

r— >0 

From Eqs. ^| and [7| it follows that F qq as well as F\ will show Coulomb-like 
behaviour at short distances and up to a normalization constant will approach the 
zero temperature heavy quark potential. We note, that this power-like short distance 
behaviour of the free energies is quite different from the behaviour of the colour 
averaged free energy at large distances (rT >> 1) where the dominant scale for the 
running coupling is set by the temperature (g(r,T) = g{T)) and high temperature 
perturbation theory is used to show that the leading order contribution to F qq is 
given by two-gluon exchange || ||. In general the statement is that 



AF qq (r,T) ( AF x (r, T) x 2 



T 



for rT >> 1 , (8) 



where ATj = Tj(r, T) — F^ with i — 1, qq. In a somewhat loose notation it thus 
often is argued that F qq /T ~ (Tx/T) 2 || |(| at large distances. In the spirit of 
our previous discussion this statement, however, has to be formulated a bit more 
carefully as Too in general will not be zero. When fixing the overall normalization 
of the free energies at short distances one no longer has the freedom to assume that 
they approach zero at large distances. In this limit the colour averaged as well as 
the singlet free energy will approach a finite value, Too(T). Using this asymptotic 
value for the heavy quark anti-quark free energy we now can define the renormalized 
Polyakov loop, 

L ren (T)=exp(-T 00 (T)/2T) . (9) 

It is obvious that Z/ en (T) is zero below T c on a lattice with infinite spatial extent, 
as T^ w iH be infinite in the confined phase. Moreover, as the renormalization 
performed at short distances only involves zero temperature physics it does not 
influence the non-analytic structure of Too which will show up in the vicinity of 
T c , the renormalized Polyakov loop. L ren thus will also reproduce all the universal 
properties extracted in the past from unrenormalized Polyakov loops 0. 



3 Numerical results for heavy quark free energies 

We have calculated the colour singlet and colour averaged free energy of a heavy 
quark anti-quark pair in a SU(3) gauge theory on lattices of size N% x N T with 
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Figure 1: The difference of the colour averaged and singlet heavy quark free energies 
calculated at two temperatures above T c calculated on lattices of size 32 3 x N T . The 
free energies have been assumed to coincide in the limit of infinite separation between 
the quark and anti-quark sources. 



N a = 32 and N T = 4, 8 and 16. For our numerical calculations we used a tree- 
level improved gauge action which includes the standard Wilson plaquette term 
and planar six link loops. This action has previously been used extensively for 
the calculation of thermodynamic quantities ||. These calculations also provided 
detailed information on the non-perturbative /3-function determined through the 
string tension aa 2 . We use this information as well as the relation T c /yfa = 0.635 
H to fix the temperature scale. In order to get access also to the singlet potential 
we have fixed the Coulomb gauge. 

Most of our calculations have been performed in the deconfined phase of the 
SU(3) gauge theory. The colour averaged and singlet free energies, calculated ac- 
cording to the definitions given in Eqs. [1] and f|, approach a constant value at infinite 
distance. As argued above the free energy should not depend on the relative orienta- 
tion of the quarks in colour space and we thus have fixed the relative normalization 
of the free energies such that they coincide in this limit. We then can analyze the 
spatial dependence of the difference (Fg q (r,T) — Fi(r,T))/T. For large distances 
this difference now vanishes by construction and we should find in the short dis- 
tance limit, rT — > 0, the relation given in Eq. [?[ As shown in Figure [I] the difference 
indeed approaches the value In 9 as expected. This also shows that the colour av- 
eraged potential at short distances indeed is dominated by the singlet contribution. 
For T > T c and distances rT < 0.1 we find that the difference deviates by less than 
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Figure 2: The shifted colour averaged (T av = F qq — Tin 9) and singlet (T x ) heavy 
quark anti-quark free energies calculated at T = 1.5 T c on lattices with temporal 
extent N T = 8 and 16. The free energies and the spatial separation r have been 
expressed in units of the square root of the string tension. The free energies have 
been matched to the zero temperature heavy quark potential (solid line) at the 
shortest distance available. The arrow points at the value (T^ — Tin 9)/ \/cr, with 
Too determined from the colour singlet potential. 



10% from the asymptotic value, In 9. We also note that the shortest distance that 
can be resolved on a lattice with finite temporal extent N T is rT = 1/N T . In order 
to analyze the short distance behaviour of the potential for rT = 0.1 it thus was 
mandatory to use the large temporal lattices which we have used here for the first 
time to study the finite temperature heavy quark free energies. 

As discussed in the previous section we expect that at short distances the colour 
averaged heavy quark anti-quark free energy is dominated by the singlet contribution 
and, moreover, approaches the zero temperature heavy quark potential, V qq {r). In 
order to compare our finite temperature calculations of the free energies with V qq (r) 
we first have to specify the latter. Using lattices with small lattice spacing V qq {r) 
has recently been calculated for the SU(3) gauge theory for distances larger than 
0.05fm and the results have been extrapolated to the continuum limit ||. For 
distances larger than a scale r , which is defined through the slope of the heavy 



quark potential |L0 



(10) 



it is known that Vq q is well described by the simple linear confining potential 
corrected by a Coulomb-like term arising from string fluctuations, i.e. V qq (r) = 
or — n/12r for r > r$ with 



ar 2 Q = 1.65- tt/ 12 . (11) 



For distances smaller than r$ we use a polynomial interpolation of the lattice data of 
Ref. || normalized such that the resulting potential smoothly joins the confinement 
potential for r > vq, i.e. we fix the free constant in the lattice results such that 
r oVqq( r = r o) = 1-65 — |. In some cases we also need the zero temperature potential 
at distances smaller than O.lro — 0.05fm. Here we use the perturbative 3-loop 



calculation of the potential in so-called qq scheme |TTJ which agrees well with lattice 



calculations up to distances 0.25r |TTJ. For the comparison of free energies with 
the zero temperature potential we give them in units of \fa which is straightforward 
since r$ and y/a are related through Eq. |TT[ This zero temperature potential is 
shown in Figure |2] as a solid line. 

For the matching of free energies to V qq we also use calculations of Fi(r,T) and 
Fq q (r,T) at off- axis separations on the lattice. At short distances these calculations 
show violations of rotational symmetry, also known from other studies of short 
distance properties on the lattice. As the matching is done at short distances we 
thus should try to adjust for these lattice artefacts. Following Ref. we replace 
Fi,qq{r) by F lm {ri) where rj 1 ^) = 4vr exp(ik-f)D$ (fc), with D$ denoting 

the tree level lattice gluon propagator evaluated in JT2]| . We thus replace the lattice 
separation r by the separation 77 which corrects for the violation of the rotational 
symmetry in the Coulomb potential calculated on the lattice. 

A comparison of the colour averaged free energy (Fq q /y/a), the singlet free energy 
{Fx/ \fo) and the heavy quark potential (Vq q / \/a) is given in Figure |2|. For clarity we 
show in this figure the shifted colour averaged free energy, F av (r, T) = Fq q (r,T) — 
Tin 9, which should coincide with the heavy quark potential at short distances. 
From Figure |2| it is evident that indeed the shifted colour averaged free energy 
and the singlet free energy approach the zero temperature heavy quark potential 
at short distances. The good agreement of F w jyfa and F\j ' \fa at short distances, 
furthermore, supports our assumption that the singlet and colour averaged {Fq q j \fo 
rather than F w j yfa) free energies approach the same constant at large distances. 
We also note that the color singlet potential calculated for N T = 8 and A^ r = 16 
agree well with each other indicating that residual cutoff effects are small. Similar 
results hold for all temperatures examined by us. 
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Figure 3: The colour singlet heavy quark anti-quark free energy in units of \/a 
versus r^fa calculated on lattices of size 32 3 x 8 for several values of the temperature 
above T c . The corresponding gauge couplings are given in Table 1. The free energies 
have been matched to the zero temperature heavy quark potential (solid line) at the 
shortest distance available. 



4 The renormalized Polyakov loop 



The analysis of the heavy quark anti-quark free energies presented in the previ- 
ous section shows that both can be renormalized by matching their short distance 
behaviour to that of the zero temperature heavy quark potential. In this way all 
divergent self-energy contributions are removed. Moreover, we have noted that the 
asymptotic large distance behaviour of F 1 (r, T) coincides with that of F qq (r, T). The 
former thus equally well determines the asymptotic value of the free energy, F^T). 
As it is obvious from Figure |2| that it is much easier to perform the matching to 
the zero temperature potential at short distances using the colour singlet free en- 
ergy F\(r,T) we have used this free energy rather than F qq (r,T) to determine the 
asymptotic value of the free energies, F^. In particular, this allows to perform the 
matching to V qq already with finite temperature free energies calculated on lattices 
with temporal extent N T = 4. Results for the colour singlet free energy calculated 
at several values of the temperature above T c are shown in Figure |3]. From these 
renormalized free energies we have determined Foo(T) at different values of the tem- 
perature. Actually, in order to avoid any fits, we use the value of F\(r, T) at maximal 
on-axis separation of the quark anti-quark sources on a lattice with spatial extent 
A^-, i.e. we define F^iT) = F 1 (N a /2,T). Systematic errors on the free energies 
have been estimated by matching them to V qq at distance rT = 1/N T as well as 
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Table 1: Change in free energy due to the presence of a heavy quark anti-quark pair 
in a thermal gluonic heat bath. The table gives results for (T)/T obtained in 
calculations on a 32 3 x 8 lattice and the screening radius defined through Eq. [12|. 
Errors on F (X) (T)/T include an estimate of the systematic errors (see text). The 
temperatures corresponding to the gauge couplings (3 are given in the first two 
columns. 

rT = y/2/N T . Some results are collected in Table |l]. 

As can be seen in Figure |3] the singlet free energy rapidly changes from the 
Coulomb like short distance behaviour to a constant value. This, of course, reflects 
the exponential screening of the two static charges. As this change is so rapid, we 
may use Foo (T)/T to define a screening radius r screen , which characterizes the onset 
of screening, through the relation 

^<jg (^screen) -^oo(-T') /1<~A 



The values for r screen /T determined in this way are also given in Table |l|. As can 
be seen r screen ~ 0.4 fm in the vicinity of T c ~ 270MeV and drops to 0.2 fm for 
T = 3 T c . We do expect that asymptotically the screening radius will drop like the 
inverse Debye mass, i.e. r screen ~ l/g(T)T. 

Using the values for ^^(T) we can now construct the renormalized Polyakov 
loop defined in Eq. |9|. Results are shown in Figure f|. The renormalized Polyakov 
loop is an order parameter for the deconfinement phase transition, which is well 
defined also in the continuum limit. Its magnitude is related to the free energy of 
a heavy quark placed in a thermal gluonic heat bath. It is obvious from Table || 
that Foo/T becomes negative at large temperatures and consequently L ren becomes 
larger than unity. We note, however, that U en is fixed only up to a multiplicative 
renormalization which results from fixing an arbitrary additive constant in the zero 
temperature heavy quark potential. 
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Figure 4: The renormalized Polyakov loop expectation value defined in Eq. [5] deter- 
mined from the asymptotic behaviour of colour singlet free energies on lattices of 
size 32 3 x JV T . 

5 Outlook 



In this paper we have defined the renormalized Polyakov loop by matching the free 
energy of a static quark anti-quark pair at short distances to the zero temperature 
heavy quark potential. We have shown that the renormalized Polyakov loop can be 
determined from the large distance behaviour of the colour averaged as well as the 
colour singlet free energy of the gg-pair. The approach has been used here to study 
the heavy quark free energy of the SU(3) gauge theory. It, however, generalizes 
without any difficulties to the case of QCD. 

In the temperature regime analyzed by us, T/T c < 6, the renormalized Polyakov 
loop is a monotonically rising function. It becomes larger than unity for T/T c ~ 
2.5. In the future it will be interesting to analyze the asymptotic behaviour of the 
Polyakov loop at even larger temperatures and determine its infinite temperature 
limit. In this limit L ven is expected to approach a constant. Figure |] suggests that 
this constant is close to unity. A more detailed analysis of the large temperature 
behaviour of Z/ en would also allow to make contact with perturbative calculations, 



which suggest that the asymptotic value may be approached from above [13j 



Finally we note that our normalization of the heavy quark anti-quark free ener- 
gies at short distances also opens the possibility for a new look at the heavy quark 
potential at finite temperature. Using the thermodynamic relations between en- 
tropy, energy and free energy, S = —dF/dT, U = —T 2 d(F/T)/dT, it is evident 
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from Figure 3 that Sq q (r,T) vanishes at short distances while it clearly is positive 
at large distances. It therefore will add a positive contribution to the total energy, 
Ug q (r,T) = Fg q (r,T) + TSq q (r,T), which becomes larger with increasing separation 
r. This shows that the change in energy due to the presence of a heavy quark anti- 
quark pair in a thermal bath is quite complex. In particular, its r-dependence is not 
only given by Fq q (r,T). This indicates that it may be misleading to use Fi(r,T) in 
potential models for the study of heavy quark bound states. 
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